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1 Scattering matrices and existence of trapped waves 

Systems with finitely many scattering channels are discussed in various branches of wave physics. One 
can mention, for instance, diffraction gratings (optics), waveguides with local perturbations (acoustics, 
microwave and quantum physics), etc. (few examples are shown in Fig^J. In all cases, far away from 
the scattering area both incident and scattered fields are sums over "modes" (specification of modes is 
given below). By another words, assuming some throughout enumeration of modes over all channels, any 
solution to the scattering problems asymptotically behaves as 
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Figure 1: Various systems with nontrivial scattering properties: "quantum wire", open resonator, diffraction 
grating. Dialing with the grating, one can reduce the problem to the strip whose width is equal to the grating's 
period (2ir in the figure) with the quasi-periodical boundary conditions @. All these systems can support trapped 
modes or guided waves. 

In what follows let us restrict ourselves either by the Helmholtz Aip + k 2 n 2 (x)ip = or by the 
Schrodinger — A-0 + V(x)ip = Ei/j equation and assume that for |x| Js i?o n(x) = 1 or V(x) = 0, 
respectively. Then modes u„(x) are solutions to the homogeneous problem for \x\ ^ Rq. For example, 
the waveguide modes are u^ = <f> n (%±.) e ±lA " 2 /\An, \i = \Jk 2 — \i n (or with E instead of fc 2 ), where 



(z,x\_), z,x± are local coordinates respectively along (outward to infinity) and across each channel, 



{/i„, 4> n (x±)} are eigenpairs to the problem in channel's cross-section z — const; the gratings modes are 

u^ = exp{±iA n z+ (n + a)y}/y/4:irA n , A„ = y k 2 — (n + a) (x± = y is one-dimensional here) satisfy the 
quasi-periodical boundary conditions 

2-rria „,| o I 2-izia 
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,=-„, d y u\ y=n = e'™d y u\ y= _„ (2) 



(here grating's period equals to 2w and a = —ks'mO, is the angle of incidence, see Fig. [JJ. It is also 
assumed in Q that each mode differ from zero only at "its" channel: by u^ we actually denote 

u t ^>u±(z,x±)r](z), (3) 

where r/(z) is a smooth cut-off function which equals to one inside the channel and zero elsewhere. 

Modes with numbers n such that A„ > propagate (oscillate) along the channel at which they are 
defined; we distinguish outgoing and incoming (respectively, u+ and u~ in our notations) modes by the 

sign of the energy flux density Im I ^u ) , where v is the outward normal to channel's cross-sections. At 
the threshold frequencies A n = 0, and the above definition of modes failed; following the paper £Q we 
define threshold (standing) modes by the following replacement: e ±lA " z / '-\/An "~> (1 T z ) /v2- 

The sum in Q) is normally taken over all propagating modes (including standing ones). Let the total 
number of such modes is 27V. Then the space of solutions with asymptotics (JIJ is TV-dimensional (see pQ 
again) and, consequently, there exists N linearly independent rows (cj~, . . . , c^, cf, . . . , c^f) , which form 
two N x N matrices C^ := ||c^„||. Both C^ are invertible (otherwise nontrivial solution to homogeneous 
problem with an asymptotics containing only outgoing or only incoming modes can exist). The unitary 
matrix § := (C~) C + is called "scattering matrix". The question of its numerical determination was 
discussed in a long series of publications (see, e.g., |2 Ej as a most general method). The approach 
suggested below can be, in particular, applied to that end as well (see Section|3J). However, our goal is a 
little bit wider. 

Let us agree that ImA n sj and include into the sum Q all modes corresponding to |ImA„| < 7 
for some positive 7. 1 Note that under our choice of the sign of ImA„ modes u~ are decay whilst u+ 
are grow as \x\ — > 00, n > N. Now the question of interest is: is there exist linear combination 
of solutions 0J) containing only decaying exponents? Existence of such a combination would 
indicate the solution to the homogeneous problem either localized at compact domain (trapped modes 
or bound states to waveguides and resonators) or propagating along a grating groove and decaying away 
from the grating (guided or surface grating waves). 

Suppose the new matrix C + is at our disposal and write it down as 

C(ll) C(12) 

C(21) C(22) 

where the block Crn\ is of size N x N and corresponds to the coefficients are in front of propagating 
modes. By straightforward verification one concludes that the condition 

detC (22) =0 (4) 

is sufficient for the existence of the aforementioned decaying solution. Consider (J2J as the equation for 
the wave frequency k (or particle energy E in quantum physics) . Thus we arrive to the existence criterion 
for localized solutions. 

It can be also shown that scattering matrix § is now the left-upper block of the product (C~) C + . 
Thus matrices C contain information about both scattering and trapping properties of the system. 

However, the numerical determination of the matrices C is not a trivial task since the asymptotics 
l|T|) contains exponents that grow and decay with different rates. The leading term dominates, and all 
coefficients cannot be found accurately. The key point of the approach suggested in the next Section 
is to avoid this difficulty. In the Section ^ we give a very brief review on previously known numerical 
approaches and provide new examples of trapped modes. 



J It is known correction term is (e tI^I) this case. 



2 Description of the approach 

The main computational idea is as follows. 

Consider the truncation of an infinite domain by removing infinite parts of all channels, starting from 
some distance R > Rq (i.e., truncating lines are Zj — R, j is the channel's number). Let ip be the solution 
to the auxiliary problem which inherits all conditions of the original problem (at an infinite domain) with 
the following artificial condition at Zj = R: 
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z i= R 



where £ > and coefficients c^ are arbitrary. 

Conditions J3J) came into scene as the differentiation of the asymptotics Q (positive £ is necessary 
to avoid possible coincidence of k 2 with the spectrum of the auxiliary problem). Let us now try to chose 
coefficients c^ in such a way that i\) approximately agrees with the right-hand part in JQ at Zj — R. 
Thus one can expect that for sufficiently large R the numbers c^ are good approximations to c„ . 

By this motivation we arrive to the condition on unknowns c^ : choose these coefficients in such a way 
that 
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It can be shown that the left-hand part in © goes to zero as R — > oo with exponential rate and 
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The proof of the similar convergence for the entries of the scattering matrix S, as well as the estimation 
of the constant 7, can be found in 0J|S]. 

From ©, one has the simple algorithm to obtain approximations c^ to c~ that are sought for. In 
fact, the functional (jSJ is quadratic with respect to qj. 
written as 

(a- B 

1 B* A+ 
where 



The Hermitian matrix of its coefficients can be 
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and the functions u„ are the solutions to the problems with artificial conditions Bj (u^ 



0. 



Note that the data to the auxiliary problems on the differences (u^ 



'-) don't contain growing as 



R — > 00 exponents as well as summation of different rate exponents. These problems have to be solved 
numerically by means of any appropriate method (e.g., finite element method). 

Finally, given the matrix (JJJ one can find its spectrum and select M smallest magnitude eigenvalues 
and corresponding eigenvectors hi,..., hu (here 2M denotes the size of the matrix) and put 



(hi,. 



I'M) 



(8) 



This finalize the description of the computation procedure. It is clear that the approach doesn't sensitive 
to the geometry of a problem. 



3 Examples of scattering data computation 

In this section we discuss the the application of the above approach to the computation of scattering 
data. 

1 Quantum control on electron stream 

Consider the two-dimensional domain D C R 2 consisting of a resonator (i.e., of a disk of radius po) 
which is connected to infinity by means of three straight channels (waveguides) of width d (see the left 
part in Fig^J); directions of axes of two channels (say 2nd and 3rd) are symmetric with respect to the 
direction of the 1st channel axis (the symmetry is not essential and is introduced to decrease the number 
of problem parameters) . We assume that the motion of an electron is bounded by the domain D and the 
wave function ip(x, y) vanishes at the boundary of D. 2 
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Figure 2: Left: Sketch of the model. Right: Transporting losses 



kiil 2 + M 2 



(solid line) and the transmission 



coefficient |sxa | (dashed line) versus energy E of an incident electron; handling potentials are fixed as V\ = 0, 
V2 — V3 — (1.57r) 2 ~ 22.2; the size of resonator po = 3d. 



The governing potential V is induced in the following way. 3 Let the resonator's walls A\, A%, A3 are 
charged by potentials V\, V2, V3, respectively; the handling is realized by variation of values Vi.,2,3. The 
whole system is shielded by three non-closed lines B±, B%, S3, each shield Bj consists of a segment of 
radius B and two rays directed along channel's walls (these shields are shown in Fig[3by dashed lines). 
Thus the potential V is the solution to Laplace equation in the unbounded domain D B (restricted by the 
shields -81,2,3) and satisfies the Dirichlet boundary conditions V\ A . — Vj, V\ B . = 0. It is known that 
such solution exponentially decays along the channels. It allows to accept, as an approximation to V, the 
solution to the problem in a finite part of D B that is located inside the disk of sufficiently large radius 8 
(with the zero boundary condition at new boundaries). 

We consider the scattering problem in the domain D under the energy range tt 2 < E < (2tt) 2 (in 
reduced units), that is, between first and second thresholds and, consequently, the size of the scattering 
matrix is 3 x 3 (one propagating mode at each channel). The examined energy range corresponds to 
realistic energy values (0.01 -j- 1) ev for d = (1 -j- 10) nm. 

Let the incident electron stream comes along the 1st channel, so the scattering data is collected in 
the first row ||sij||, j — 1,2,3, of the scattering matrix S. Now the question is if one can choose any 
combination of the energy E and handling potentials Vi, V2, V3 such that the transport of an electron 
stream either to 2nd or to 3rd channel is close to certain (i.e., either | S12 1 ~ 1 or |si3 j ~ 1)? If yes, the 
transport can be switch from second to third channel by trading places of V\ and V3. 

2 The Schrodinger equation is written in reduced units when unit of length is d and the unit of energy is ft? / (2m* d?) , 
m* is an effective electron mass. 

3 This model was suggested by Prof. L.M.Baskin (St. Petersburg, Russia); more detailed description of the model and 
extended results are in press [H]. 



Some of the results are given in the right-part in Fig|2J It is seen that transporting losses can be 
reduced below 0.1% 4 by variation of E (the same effect can be achieved by variation of V as well, not 
shown here). 

2 Conductance of bent waveguides 

Two-dimensional bent waveguides are often considered as models of microwave devices or quantum wires 
(see |S]). Scattering problem for such a model is to determine the matrix § corresponding to the solution of 
the free Schrodinger (Helmholtz) equation in a curved strip with asymptotics (JTJ and satisfying Dirichlct 
boundary conditions. Normally the implementation of a device presumes the strip of piecewise constant 
width (e.g., a bent waveguide is connected to open leads, as shown in the left part of FigOU; this case the 
definition of modes and scattering matrix is related to leads. 

Typically the wave process in a system is energized by an incoming single mode (number m) of inlet 
lead, and the observed quantity is the normalized conductance G m = J2 n \ s mn\ , where sum is taken over 
all outgoing modes of outlet lead. 
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Figure 3: Left: A bent waveguide (numerous dimensions of the bend are not specified) of width d connected 
with two open leads of width H; Right: Normalized conductance Gi (due to the lowest mode in the lead) versus 
dimensionless frequency kd: solid line — H/d = 2, dashed line — H/d = 4. Note that the location of conductance 
maxima doesn't depend on the ratio H/d and the first peak is observed (below the bent waveguide threshold n) 
at kd « 3.097 (compare with the trapped mode frequency for the waveguide with the same bent part, see Figlol, 

It is clear that modes can propagate throughout the leads at frequencies below the threshold (kd = tt) 
of the bent strip. Transmission of energy at such frequencies through the bend is classically forbidden. 
However, that becomes possible due to quantum effects. The computation results are shown in the 
right part of FigO the first conductance peak is observed below the threshold; above the threshold, the 
conductance exhibits an oscillating behavior. The frequency corresponding to the first conductance peak 
is close to the trapped frequency of the infinite bent waveguide (see Section 0J and can be thought as 
resonance frequency. 



4 Computation of trapped modes: brief review and new exam- 
ples 

The localized solutions have been attracted a lot of attention during last decades as examples of nonuniquc- 
ness to a scattering problem which, at the same time, exhibit abnormal physical properties. Depending 
on the area of application these solutions are either called trapped modes (microwave and/or acoustic 
waveguides, water waves), or bound states (quantum wires and photonic crystals), or guided (Rayleigh- 
Bloch) waves (diffraction gratings). The framework of a workshop publication doesn't suit for detailed 

4 The accuracy of computation of the scattering coefficients better than 0.005; thus, the total transporting losses, i.e., 
| si] 



I-S13 1 2 , exceeding (2 



-5)10 5 were found safely. 



review of tens of papers devoted to the subject. The whole spectrum of papers can be conventionally 
divided into three categories: proof of existence or non-existence of trapped modes, asymptotic estimate 
of trapped eigenfrcquencies (energies) and numerical approaches. Not mention the first category com- 
pletely, we point out the key names and selected papers related to the last two categories. The more 
comprehensive review in a part related to the quantum applications can be found in the book \7\. 

Concerning waveguide problems, the asymptotic results can normally be obtained if the geometry 
is a perturbation of a straight strip. This may be a case for small indentation [H], laterally coupled 
waveguides through small window(s) [§J, slightly curved strips and tubes |l()j . etc. (see, e.g., |llj for 
asymptotic results in gratings). Normally, the asymptotic formulas match numerical computation for 
some range of "small" parameter; at the same time, numerics is often fails when perturbation becomes 
very small. 

As far as known to the author, the overwhelming majority of papers dedicated to the numerical 
detection of trapped modes dealt with objects of a relatively simple geometry. The typical methods 
were either the separation of variables in sub-domains with subsequent matching of infinite series, or the 
reduction of a problem to the integral equation (the explicit knowledge of Green's function is required), or 
the application of a variational technique (the sufficiently reach set of trial functions satisfying boundary 
conditions is required). All these methods were successfully utilized and a wide variety of trapped modes 
found by research groups represented by papers E3 [El EI! ■ 

However, all the mentioned methods are of limited universality and applicability range. The developed 
approach allows us to contribute into the scope providing successfully found new examples of trapped 
modes. Due to the generality of the approach the list of examples can be very wide, below are few of 
them. 5 

The Dirichlet trapped modes eigenfrequencies for the straight waveguide with a local indentation are 
shown in Fig0J Such model was the subject of rigorous studies in [H] where existence of trapped modes 
below the threshold were proved and asymptotic estimate of them found. In Fig^Jthe numerical results 
are given in comparison with the asymptotic formulas of 8 . It is seen that both results match each other 
although the range of applicability of asymptotic formulas is relatively narrow. 





* 


fc> 




. 










d 




R o 


■ 












Figure 4: Left: Waveguide with an one-side indentation. Right: The Dirichlet mode trapped eigenfrequencies kd 
(solid line) versus the length b/d of the indentation (whilst Ro/d = 1) in comparison with the asymptotic results 
of the paper [H] (dashed line); the inset shows the vicinity of the threshold in an expanded range. 

Examples in Fig[5] concern the trapped modes for that no numerical estimates were known. In the 
left part we show the map of \iptrappedix, y)\ for a bent waveguide which bent structure is the same as in 
the example of FigEl (but with no leads attached). After replacing of the waveguide's straight parts by 
wider leads the eigenfrequency (energy) of this trapped mode moves to a resonance in the complex plane. 
This resonance works as a "bridge" providing the sub-threshold conductance peak in FigEI 

In the right part of FigJSJwe provide an example of trapping properties of the domain with relatively 
complicated geometry. The domain in question consists of three smoothly and symmetrically connected 
channels (the symmetry and the number of channels are not essential); some disk is taken out around 



5 The Workshop presentation contains more ones, including examples related to the diffraction gratings. 
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Figure 5: Left: The intensity map of the trapped mode in a bent waveguide; the bent part is the same as in 
FigEl trapped mode frequency is kd fa 3.091 (which is slightly below the sub-threshold conductance frequency 
for the bent waveguide with adjoint leads, see Section [5J. Right: The dimensionless trapped eigenfrequencies kd 
(versus the size of disk hole) for the domain shown in the inset. 

the symmetry center. The Dirichlet boundary conditions are assumed at the outdoor domain boundaries 
whilst the Neumann conditions are set at the interior circle boundary. This domain with no hole (a/d = 0) 
can be thought as a bent waveguide with an "infinite" indentation, thus the existence of the Dirichlet 
trapped mode below the threshold can be clearly predicted, and it is numerically found as kd « 2.594. 
The trapped eigenfrequencies approaches the threshold versus increasing the disk hole radius. 

5 Conclusion 

The general approach for numerical computation of scattering and trapping properties of systems with 
finitely many channels to infinity quantum, microwave and optics physics was discussed. The approach 
can be adopted for the numerical search of resonances for the same type of scattering systems; this work 
is under development. 
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